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Résumé : 

De nombreux instruments de musique (comme les instruments à vent ou à cordes frottées) ont un fonctionnement 
reposant sur l'entretien d' oscillations par une action continue du musicien. Dans le cas de la clarinette, l'anche, 
vibrant au gré des ondes acoustiques dans le corps de l'instrument, module l'apport d'énergie qu 'est le souffle du 
musicien. Nous nous intéressons au seuil d' oscillation, valeur de la pression dans la bouche à partir de laquelle 
^ ■ un son peut prendre naissance et cherchons à comprendre comment le comportement dynamique de l'anche, par 

\ sa collaboration avec les résonances du tuyau, influe sur ce seuil. Nous présentons ici une méthode reposant sur 



o 

m 



la décomposition modale du résonateur acoustique qui autorise une étude au-delà du seuil, et la confrontons avec 
une résolution directe du problème de recherche de seuil. 

Abstract : 

Many musical instruments, as for example woodwind instruments, flûte or violins, are self-sustained oscillating 
Systems, i.e. musician enacts as a continuous energy source to drive an oscillation in the passive resonator, 
the body of the instrument, by means of a nonlinear coupling. For single reed instruments like clarinet, there 
exists a minimal value of mouth pressure beyond which sound can appear. This paper deals with the analysis 
. of this oscillation threshold, calculated using a modal décomposition of the resonator, in order to have a better 

' compréhension ofhow reed characteristics, such as its strength and its damping, may influence the attack transient 

\^ , of notes. 
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^ ■ acoustique musicale, décomposition modale, seuil d'oscillation 



1 Introduction 

Dans les instruments de musique à anche simple tels que la clarinette comme dans d'autres 
systèmes comprenant un guide d'onde muni d'une valve libre à une de ses extrémités, des oscil- 
lations peuvent apparaître dès lors que la pression d'alimentation dépasse une valeur critique. 
Ce seuil d'oscillation, directement lié à la notion de facilité d'émission du son pour les ins- 
truments de mus i que, a fait l'objet d'un grand no mbre d'études, aussi bien expérimentales (cf . 
Dalmont et coll.\ (l2006|) . IWilson et Beavers qu'analytiques (cf. I Palmont et colÛ |2005|) . 



Grand ( 19941) ) sur des modèles simplifiés, comme par exemple le modèle dit de Raman. La 



motivation principale de ce travail est de proposer une méthode numérique d'investigation des 
seuils d'oscillation reposant sur l'étude linéaire de stabilité. L'originalité de l'approche consiste 
en une description modale du résonateur acoustique, généralisant ainsi des travaux précédents 
portant sur l'interaction entre un mode d'anche lippale et un mode de tuyau dans le cas des 



cuivres, cf. ICuUen et coll.] (|2000t) . l'anche étant ici décrite par son p remier mode de flexion. La 



méthode de décomposition modale a été précédemment utilisée par DebutI (|2004î) pour le calcul 



des auto-oscillations dans le domaine temporel. L'intérêt de la méthode est de pouvoir aisément 
être transposée aux différents instruments à anche, aussi bien pour les instruments coniques que 
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pour les anches doubles et les cuivres. 
2 Le modèle linéarisé 

Nous modélisons la clarinette comme le couplage d'un résonateur acoustique (le corps de 
l'instrument) et d'un élément vibrant (l'anche) par l'écoulement d'air à travers le canal d'anche 
sous l'effet de la différence de pression entre la bouche du musicien et l'entrée du tuyau. 

Un des objectifs du travail est de caractériser la manière dont le comportement dynamique 
de l'anche intervient dans le seuil d'oscillation du système couplé. Nous limitons ici cette dy- 
namique au premier mode de flexion de l'anche, le mouvement de l'anche pouvant ainsi être 
décrit par une équation différentielle du 2°ordre liant le déplacement de son extrémité libre y{t) 
(par rapport à sa position au repos yo) à la différence entre la pression dans la bouche pm et dans 
le bec p(t) : 



UJ. 



où K est un coefficient surfacique de raideur qui peut être évalué à partir d'un essai statique, cUr 
et Qr étant la pu lsation propre et le coeffi cient d'amortissement de l'anche estimés à partir d'un 



essai de lâché ( Wilson et Beaversl (|1974t) ). 

Le couplage entre vibration de l'anche et ondes acoustiques dans le corps de l'instrument 
est réalisé de manière aéro-acoustique par l'écoulement de l'air de la bouche vers le bec de la 
clarinette à traver s 1 ' ouverture entre l'anche et la table. L'utilisation du théorème de BemouUi 
a été justifiée par Hirschber^ ( 1995 ) sous les hypothèses d'écoulement stationnaire et incom- 



pressible (masse volumique p), de grande différence de section entre la bouche et le canal et 
de dissipation de l'énergie cinétique du jet par turbulence à la sortie du canal. En supposant de 
plus l'ouverture rectangulaire de largeur W et de hauteur y{t), il permet d'établir la relation 
non-linéaire entre le débit u{t) entrant dans l'instrument et la différence de pression : 



u 



(t) = Wy{t) ^sign {pm - p{t))^/\pm - p{t) \ si y(t) > 0, sinon «(t) = 0. (2) 



L'étude se portant au voisinage du seuil d'oscillation afin de déterminer quels sont les modes du 
système couplé qui naissent en premier, nous nous plaçons dans le cadre d'une étude linéaire 
de stabilité, ce qui autorise la linéarisation de la caractéristique de couplage et le passage dans 
le domaine fréquentiel (avec une convention de dépendance temporelle en expjcjt) pour de 
faibles écarts au régime statique. Il vient alors : 

U{u) = ^ W(u;) _ ^ A I^l") p(^) et Y{uj) = ^ \ PM (3) 



L'admittance aéro-acoustique (réduite) vue par le résonateur acoustique est donc : 

où Zc = pc/S est l'impédance caractéristique du tuyau (liée à sa section S" et à la célérité c du 
son en espace libre), 7 et ^ sont deux paramètres adimensionnés représentant la pression dans 
la bouche et la « pince »du musicien : 



7 = 1^ et C = ^W% (5) 
Kyo V pK 
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Le résonateur acoustique est décrit par son impédance d'entrée Zf., en faisant l'hypothèse 
d'un tuyau cylindrique de longueur L et de rayon r ouvert à son extrémité (l'impédance de 
rayonnement est considérée nulle), et en prenant en compte un amortissement des ondes a par 
des phénomènes visco-thermiques aux parois. Cette impédance peut être décomposée sur les 

premiers modes propres du résonateur passif définis par les pulsations propres Vtn (ce sont 
celles du tuyau sans pertes légèrement abaissées du fait de la dispersion résultant des pertes 
visco-thermiques) et leurs coefficients d'amortissement a„ (estimées à partir de la valeur de a 
à la pulsation propre associée) : 



c étant la célérité du son en espace libre. 

3 Modes propres du système linéarisé et seuils d'oscillation 

L'équation caractéristique du système linéarisé est obtenue en égalant impédance aéro- 
acoustique et impédance d'entrée du résonateur. Avec la formulation modale de l'impédance 
du tuyau, il vient : 



La démarche proposée est celle de la recherche du seuil d'oscillation par les méthodes 
d'étude de stabilité linéaire. La recherche des solutions joj de l'équation caractéristique pour 
une configuration donnée (L, r, ujr, Qr, 7 et ^ fixés) permet d'obtenir les modes propres du 
système linéarisé autour du régime statique. Les parties imaginaire et réelle correspondent res- 
pectivement à la fréquence et à l'amortissement des différents modes propres : 



Un des avantages de la formulation modale de l'impédance d'entrée du résonateur est de 
pouvoir faire apparaître l'équation caractéristique sous forme d'un polynôme en ju, ce qui per- 
met l'utilisation d'algorithmes efficaces pour obtenir les (2N-\-2) solutions réelles ou complexes 
conjuguées. La question du nombre de modes considérés, c'est-à-dire de l'ordre de troncature 
du développement modal de la formulation est examinée plus bas. 

Il est possible d'observer l'évolution des fréquences propres complexes en fonction de 7, 
image de la pression dans la bouche, la recherche des modes propres étant répétée pour chaque 
valeur de 7. Sur la figure [H on peut visualiser une première plage de pression où toutes les 
solutions sont amorties ( 3îe {{juj)) < 0) : le régime statique est le seul régime stable pour de 
faible pression d'alimentation. Il existe ensuite une valeur de 7 pour laquelle une unique solu- 
tion a une partie réelle nulle, les autres étant à partie réelle négative : c'est le seuil d'instabilité 
du régime statique. Pour des valeurs supérieures, il existe au moins un mode propre amplifié : 
le régime statique est instable. Dans ce dernier cas, les oscillations peuvent alors croître expo- 
nentiellement jusqu'à ce que la non-linéarité contenue dans la caractéristique de débit vienne 
saturer l'amplitude des oscillations. 

Sur la figure [Hest représentée l'évolution des modes propres dans le cas où cUrL/c = 3. Les 
fréquences propres du système couplé ont des valeurs proches des fréquences de résonance du 
tuyau isolé (a; ~ (n — l/2)nc/L) et de l'anche seule (tu = 1), les fréquences propres variant 



2c 



N 





(7) 



(8) 



18 Congrès Français de Mécanique 



Grenoble, 27-31 août 2007 




faiblement avec la pression d'alimentation. Le seuil d'instabilité du régime statique est de 7 ~ 
0.3, avec une oscillation émergente à une fréquence correspondant à la première résonance du 
tuyau. Dans le second cas (figure l2l uJrL/c = 6), le seuil, obtenu pour une valeur de pression 
proche de la précédente, est associé à la seconde résonance du tube, celle-ci correspondant à une 
fréquence légèrement inférieure à la résonance d'anche. L'hypothèse d'un comportement de 
l'anche réduit à sa raideur aboutirait à un seuil associé au premier mode du tuyau inférieur aux 
seuils des modes supérieurs du fait de l'augmentation des pertes avec la fréquence, cependant 
l'interaction entre les résonances acoustiques avec la résonance mécanique est responsable de 



cette i nversion de seuils entre premier et second mode, comme expliqué par IWilson et Beavers 
(Il974h . 

L'ordre de troncature de la décomposition modale a été choisi selon un critère lié à la réso- 
lution de notre algorithme de recherche de seuil et à la vitesse de convergence du seuil 7 avec 
le nombre de modes pris en compte. La troncature à = 20 modes a été retenue pour les 
simulations numériques, elle correspond à une modification du résonateur engendrant une dé- 
viation du seuil d'oscillation de moins de 0.01, qui est la tolérance utilisée lors de la recherche 
de zéros. 



Comparaison avec une résolution directe 

La détermination du seuil d'oscillation avec l'approche modale a été réalisée numérique- 
ment avec la méthode de la position fausse en s'intéressant à la fonction : 

F : 7 — > max { 3îe (juj) / jco solution de O}, (9) 

le changement de signe de cette fonction correspondant à la perte de stabilité du régime statique. 

Nous avons confronté (cf. Table [B les résultats obtenus par cette méthode avec ceux issus 
d'une résolution directe de l'équation caractéristique en (9, 7) avec 9 = cu/cUr à valeur réelle. 
Une des difficultés de la résolution directe est liée à la nature transcendentale de l'équation 
caractéristique avec la formulation de l'impédance d'entrée en tan . Le nombre de solutions 
étant infini, il faut s'assurer que la recherche des solutions amène bien au seuil le plus bas. 
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FiG. 2 - Évolution des fréquences propres complexes en fonction de la pression dans la bouche : 

L = 32cm, r = 7mm, Ur = 2it x lOOOrad/s, qr = 0.3, C = 0.2. 
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Tab. 1 - Comparaison des seuils obtenus par l'approche modale (notée DM) et par résolution directe 
(RD) pour r = 7mm, w,. = 27r x 750rad/s, Qr = 0.4, C = 0.13. 



4 Couplage d'un mode d'anche et d'un mode du tuyau 

Dans ce cas simple, on obtient des résultats analytiques sur la manière dont l'écoulement va 
modifier les modes des systèmes. Avec l'hypothèse d'un tuyau à un seul mode et en ne tenant 
compte que des fréquences positives, une approximation de D et Zi est donnée par : 

2 L 

D^^{uj) c:^ iuj — ujr — jqriOr/2) et ZrVc;) ~ j — (a; — fi„ — ja„c), (10) 

UJr C 

ce qui, à partir de l'équation (|7]), aboutit à la relation de couplage : 

{uj-uj^){uj-uj^) ~ j^Cy? avec = Ur (l et ^+ = + 3anC + 

(11) 

On retrouve le fait que les deux modes sont couplés par l'écoulement, mais le mode acoustique 
à considérer n'est pas celui du tuyau « ouvert-fermé » : puisque 7 < 1 (sinon l'anche plaque en 
régime statique), les pertes sont augmentées, ce qui s'explique par le comportement acoustique 
résistif de l'écoulement à l'entrée du tuyau. 

Un résultat remarquable que l'on obtient sur ce modèle simplifié est que lorqu'une solution 
devient instable (l'exponentielle temporelle associée devient croissante), l'autre solution est 
toujours stable. En effet, la somme des deux solutions vaut {uj^ + cj^), et cette quantité est 
toujours à partie imaginaire positive. Dans cette configuration à deux modes, les seuils étudiés 
sont donc bien des seuils d'instabilité du régime statique. 
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5 Conclusions 

La comparaison avec une méthode de détermination directe des seuils d'instabilité du ré- 
gime statique permet de valider la démarche proposée dans cette communication. Il est possible 
d'expliquer les faibles écarts constatés pour 7 par le caractère itératif de la méthode utilisée, 
avec une condition d'arrêt portant sur un encadrement de la solution à ^7 = 0.01. Un autre 
facteur est la troncature de la décomposition modale de l'impédance d'entrée du résonateur. En 
effet, la troncature correspond à une modification du guide d'onde, d'autant moins significative 
que le nombre de modes pris en compte est grand. Toutefois, l'écart entre les résultats des deux 
méthodes est suffisamment faible pour ne pas chercher à travailler avec un plus grand nombre 
de modes. 

Un des intérêts majeurs de cette démarche est de pouvoir être appliquée pour tout résonateur 
dont l'impédance d'entrée peut se mettre sous forme modale, puisqu'il suffit de connaître l'am- 
plitude, la pulsation et le facteur de qualité de chacun des modes, ce qui peut fait par ajustement 
de paramètres sur une mesure fréquentielle d'impédance d'entrée. Il est également possible 
d'étendre à des éléments mécaniques à plusieurs modes de flexion et/ou de torsion, à la prise en 
compte de phénomènes supplémentaires comme le débit d'anche. 

Outre l'étude des seuils d'instabilité du régime statique, l'approche proposée permet éga- 
lement de s'intéresser à ce qu'il se passe au-delà du seuil. En effet, dans la limite de validité 
de la linéarisation de la caractéristique de débit, le modèle linéarisé reste adapté à l'étude des 
démarrages des auto-oscillations, notamment en ce qui concerne la vitesse des transitoires, tout 
au moins jusqu'à ce que les non-linéarités viennent saturer l'amplitude des ondes acoustiques. 

Les travaux présentés dans cette communication ont été réalisés dans le cadre de CONSONNES, 
projet financé par l'Agence Nationale de la Recherche. 
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